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FOREWORD 


This  Final  Technical  Report  presents  the  results  obtained  under  Contract 
MDA903-79-C-0186,  DARPA  Order  Number  3700-9Y1 0-6271 2E.  The  views  and 
conclusions  contained  in  this  document  are  those  of  the  authors  and  should  not 
be  interpreted  as  necessarily  representing  the  official  policies,  either 
express  or  implied,  of  the  Defense  Advanced  Research  Projects  Agency  or  the 
United  States  Government. 

This  contract  with  the  Defense  Supply  Service  -  Washington  was  initiated  by 
the  Defense  Advanced  Research  Projects  Agency  and  was  under  the  technical 
direction  of  Dr.  Michael  J.  Buckley,  Program  Manager,  Material  Sciences 
Office,  DARPA,  and  Dr.  Henry  Graham,  Air  Force  Materials  Laboratory/LIA!. 

The  Program  Manager  for  Adaptronics  was  Dr.  Anthony  N.  Mucciardi,  the 
Principal  Investigator  was  Mr.  Dixon  Cleveland,  and  major  contributors  to  the 
work  were  Messrs.  Peter  M.  Garafola  and  Basil  A.  Decina. 

The  reaction-bonded  silicon  nitride  process  data  used  in  this  work  was 
obtained  by  and  provided  to  Adaptronics  by  the  AiResearch  Manufacturing 
Company  and  AiResearch  Casting  Company  which  are  Divisions  of  the  Garrett 
Corporation.  Garrett  obtained  the  data  from  experimental  work  performed  under 
Phase  2  of  an  Air  Force  Materials  Laboratory  project  [1]  to  demonstrate 
capability  of  increased  yield  of  slip-cast  ceramic  vanes  as  components  for 
high-performance  turbine  engines.  Mr.  David  W.  Richerson  of  AMC  and  Michael 
E.  Rorabaugh  of  ACC  were  primarily  responsible  for  the  data  collection  and 
compilation. 


ii 


A  program  was  conducted  to  model  the  modulus  of  rupture  (MOR)  strength  using 
Adaptive  Learning  Networks  (ALN's)  for  aircraft  engine  components  oroduced  by 
a  slip-cast,  reaction-bonded,  silicon-nitride  production  process.  The  primary 
objectives  of  the  'work  were  to  identify  key  process  variables  and  to  predict 
optimum  values  for  those  variables  as  a  guide  for  further  experimentation. 
Nonlinear  models  have  been  synthesized  that  predict  MOR  with  an  average  error 
of  about  4  ksi  over  a  range  from  18.6  to  47.8. 


The  manufacturing  and  analysis  work  done  to  date  has  demonstated  the  feasi¬ 
bility  of  modeling  the  slip-cast  RBSN  process  with  the  Adaptive  Learning 
Network  methodology  and  is  viewed  as  the  first  iteration  in  the  optimization 


procedure  which  is  ultimately  aimed  at  finding  those  manufacturing  conditions 
which  will  produce  the  strongest,  most  consistent  material  strengths. 
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1 .  INTRODUCTION 


This  work  involves  the  Adaptive  Learning  Network  (ALN)  modeling  of  a  slip- 
cast,  reaction-bonded  silicon-nitride  (RBSN),  ceramic  production  process.  The 
primary  objectives  of  the  work  are  to  identify  the  key  parameters  of  the 
production  process,  to  seek  an  optimum  set  of  parameter  values,  and  to  develop 
control  procedures  which  yield  highest  and  most  consistent  part  strengths. 

1.1  THE  ROLE  OF  MODELING  IN  PROCESS  OPTIMIZATION 

Optimization  of  the  ceramics  manufacturing  process  involves  an  experimental 
search  for  the  ideal  manufacturing  conditions.  Typically  an  iterative 
procedure  is  employed.  At  each  round  or  iteration  of  the  search  procedure, 
several  manufacturing  experiments  are  performed.  The  experimenter  then 
analyzes  the  results  and,  based  on  this  analysis,  formulates  hypotheses  on 
parameter  regions  which  are  projected  to  yield  improved  performance.  These 
projections  form  the  basis  for  designing  the  next  round  of  experiments.  The 
experimental  results  at  each  round  serve  to  confirm  or  reject  the  hypotheses 
formed  in  earlier  rounds  and  will  influence  further  experimentation 
accordingly. 

The  efficiency  and  effectiveness  of  the  search  procedure  is  highly  dependent 
upon  the  adequacy  of  the  analysis  function.  Implicit  in  any  analysis  is  the 
formulation  of  a  model  of  the  process.  The  model  is  the  tool  which  gives  the 
ability  to  predict  future  performance  of  the  process  based  upon  past 
observations. 

Modeling  has  three  primary  roles  in  process  optimization.  First,  it  is  used 
to  guide  the  experimental  search  such  that  a  minimum  amount  of  experimental 
time,  cost,  and  effort  are  incurred  in  converging  on  the  optimum.  Second,  the 
modeling  will  reveal  which  variables  are  important,  what  their  optimum  values 
are,  and  to  what  tolerances  they  must  be  controlled.  Third,  if  certain 
important  variables  have  not  been  measured,  model  performance  metrics  will 
indicate  that  crucial  information  is  missing,  which  will  guide  the 
experimenters  in  identifying  additional  variables  to  be  instrumented  in 
subsequent  experiments. 


1.1.1  Guiding  the  Search 

Due  to  the  expense  and  time  of  performing  the  laboratory  experiments,  it  is 
desired  that  each  experiment  provide  a  maximum  amount  of  new  information  about 
better  manufacturing  conditions.  The  procedure  for  designing  the  next 
experiment  consists  of  two  steps.  First,  the  model  is  updated,  or  retrained, 
using  all  the  data  gathered  to  date.  Secondly,  the  model  is  interrogated,  or 
searched,  for  regions  of  highest  predicted  strength.  The  search  is  permitted 
to  range  somewhat  beyond  the  region  of  data  already  collected,  and  typically 
the  highest  areas  are  outside  the  data  regions.  Though  the  model  is  not 
expected  to  be  highly  accurate  outside  the  data  regions,  it  is  a  good 
hypothesis  that  collection  of  further  experimental  data  in  the  high  areas  will 
yield  the  most  useful  results.  Subsequent  experiments  are  then  performed  in 
these  projected  high  regions. 


The  optimization  process  is  complete  when  the  peak  operating  conditions  have 
been  found.  Three  conditions  must  be  satisfied  to  verify  the  peak.  First, 
the  model  surface  must  exhibit  a  major  peak  or  high  region  where  any 
deviations  from  that  area  result  in  lower  predicted  strength.  Secondly,  all 
the  regions  on  and  surrounding  the  peak  should  be  well  supported  by 
experimental  data.  Third,  the  accuracy  of  the  model,  as  judged  on  the  data 
immediately  surrounding  the  peak,  should  be  quite  high,  indicating  that  the 
model  properly  accounts  for  all  key  variables  and  adequately  represents  the 
process. 


1.1.2  Determinincr  Controls  and  Tolerances 


Once  the  optimization  is  complete,  the  optimum  values  of  the  manufacturing 
parameters  may  be  found  by  locating  the  model  peak.  Control  tolerances  are 
then  determined  by  varying  each  parameter  from  its  optimum  value  and  observing 
the  consequent  effect  on  predicted  strength. 

Nonlinear  variable  interaction  in  the  peak  region  should  also  be  examined  to 
determine  whether  a  deviation  on  one  parameter  significantly  varies  the 
tolerance  on  another.  Trade-off  analyses  may  be  undertaken  _.o  determine 
whether  some  non-optimum  set  point  would  produce  less  parameter  sensitivity 
and  thus  more  reliable  results.  These  control  analyses  are  performed  by 
investigation  of  models  and  do  not  require  further  laboratory  experimentation. 


1.1.3  Finding  Missing  Information 

If  one  or  more  crucial  process  parameters  are  not  measured  in  the  production 
experiments,  there  is  no  way  for  the  modeling  algorithm  to  identify  them 
specifically.  But,  by  an  analysis  of  the  accuracy  of  the  models,  it  can  be 
determined  whether  further  information  is  needed  to  model  the  process 
accurately,  and  an  estimate  of  the  importance  of  that  missing  information  can 
be  obtained  from  the  model  errors.  Knowledge  of  the  need  for  further 
information  is  ver”  useful  to  the  experiment  designer  as  he  seeks  to  identify 
additional  variables  which  must  be  measured.  Highly  accurate  model 
performance  is  an  indication  that  all  the  meaningful  variables  have  been 
accounted  for. 


1.2  REQUIREMENTS  OF  THE  MODELING  ALGORITHM 

The  effectiveness  of  a  process  optimization  procedure  is  highly  dependent  upon 
the  power  of  the  modeling  algorithm.  There  are  four  key  requirements  on  the 
modeling  algorithm  which  are  generally  not  fully  met  by  conventional  modeling 
approaches  but  which  are  met  by  the  Adaptive  Learning  Network  methodology. 


1.2.1  Nonlinearity 


When  modeling  a  process  for  the  purpose  of  optimization,  the  final  model  must 
embody  a  convex  surface  to  represent  the  region  of  the  optimum  peak.  A  linear 
representation  is  not  adequate  because  it  has  no  finite  optimum.  In  a  poly¬ 
nomial  expansion,  a  convex  surface  can  be  no  less  than  a  quadratic  in  each 
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input  variable,  and  in  many  instances  the  degree  may  be  higher.  Processes 
such  as  ceramics  manufacturing  are  very  likely  to  have  subtle  but  distinct 
nonlinear  interactions  between  variables.  The  ALN  method  automatically 
considers  higher  order  terms  and  nonlinear  interactive  terms  for  all  candidate 
input  variables. 


1.2.2  Automatic  Synthesis  of  Model  Structure 

In  most  conventional  modeling  approaches,  the  user  selects  the  mathematical 
structure  of  the  model,  or  a  small  set  of  possible  structures,  and  the 
algorithm  determines  the  coefficient  values  which  produce  the  best  fit  to  the 
data.  Specifying  a  model  structure  for  a  process  as  complex  as  a  ceramics 
manufacturing  process  is  very  difficult,  and  if  the  proper  model  structure  is 
not  chosen,  the  model  accuracy  will  be  poor  no  matter  how  good  or  complete  the 
data  is. 

The  ALN  model  synthesis  algorithm  automatically  generates  the  model  structure 
as  well  as  the  coefficient  values.-!/  The  routine  generates  its  models  by 
systematically  incorporating  only  those  polynomial  terms  and  functions  which 
provide  the  maximum  performance  improvement  with  minimum  increased  model 
complexity.  The  structure  synthesis  procedure  thus  automatically  selects  the 
most  important  process  variables  and  defines,  in  mathematical  terms,  their 
relationships  to  material  strength. 


1.2.3  Prevention  of  Overfit 


A  trained  model  is  said  to  be  "overfitted"  if  it  produces  small  errors  on  the 
data  upon  which  it  was  trained  but  performs  poorly  on  similar  data  that  were 
not  used  in  training.  An  overfitted  model  is  thus  not  useful  as  a  predictive 
tool  which  can  forecast,  with  acceptable  accuracy,  what  the  results  of  a 
future  experiment  will  be. 

In  process  optimization,  where  prediction  of  :■  results  of  future  experiments 
is  important  to  the  minimization  of  the  number  of  experiments  which  are  to  be 
performed,  and  in  control  synthesis,  where  variable  sensitivities  must  be 
accurately  estimated,  overfit  must  be  minimized.  Overfit  generally  occurs 
when  the  model  is  more  complex  than  is  statistically  justified  by  the  given 
data  base. 

On  the  other  hand,  it  is  desirable  to  obtain  as  much  information  from  the  data 
as  possible  to  support  a  model  of  a  complex  process.  The  ALN  algorithm 
employs  information  theoretic  measures  which  permit  the  growth  of  model 
complexity  up  to  but  not  beyond  the  extent  justified  by  the  given  data  base. 


— ^The  ALN  structure  uses  the  form  of  a  polynomial,  rather  than  an  exponential 
or  some  other  transcendental  form?  however,  the  polynomial  expansion  is  a 
very  powerful,  general  representation  that  can  mathematically  approximate 
any  continuous  function. 


1.2.4  Treatment  of  Limited  Data  Bases 


Most  conventional  modeling  approaches  require  more  observations  than  there  are 
variables.  In  the  early  stages  of  ceramics  optimization,  this  situation  does 
not  exist.  It  is  desirable  to  vary  on  the  order  of  100  factors  but  there  may 
be  only  25  to  50  experimental  observations.  In  conventional  experimental 
designs,  only  a  few  parameters  are  varied  while  all  others  are  held  con¬ 
stant.  The  ALN  method  can  extract  meaningful  models  from  a  limited  data  base, 
even  when  the  number  of  varying  parameters  far  exceeds  the  number  of  experi¬ 
mental  observations. 

In  summary,  the  ALN  modeling  procedure  is  specifically  suited  for  the  non¬ 
linearity,  unknown-structure,  non-overfit,  and  limited-data  requirements  for 
optimization  of  complex  processes. 

1.3  SUMMARY  OF  WORK  PERFORMED  TO  DATE 

Laboratory  work  was  performed  by  the  Garrett  Corporation  on  a  slip-cast  RBSN 
process  to  produce  test  bars  under  35  different  sets  of  manufacturing  condi¬ 
tions.  For  each  condition,  the  manufacturing  parameters  were  recorded  and 
test  bars  were  destructively  tested  to  obtain  data  on  the  resulting  material 
strengths.  The  strength  parameters  were  room-temperature  modulus-of-rupture 
(MOR)  and  Weibul  modulus.  Adaptronics,  Inc.  performed  the  modeling  analysis 
of  the  data. 

In  the  course  of  the  modeling,  three  types  of  models  were  synthesized.  First, 
strength  and  strength  variance  were  modeled  as  a  function  of  the  independent 
input  variables,  such  as  slip  proportions  and  sintering  temperatures.  Second 
strength  and  strength  variance  were  modeled  as  a  function  of  the  intermediate 
process  variables,  such  as  nitrided  density  and  weight  gain.  Third,  the 
intermediate  process  variables  were  modeled  as  a  function  of  the  independent 
inputs.  The  combination  of  these  three  sets  of  models  shows  the  overall  flow 
of  effects  through  the  production  process.  Actual  material  strengths  varied 
from  approximately  19  to  48  ksi,  and  the  models  predicted  these  strengths  with 
an  average  error  of  4  ksi  over  the  total  range  of  29  ksi.  Key  process 
parameters  and  the  means  by  which  they  influenced  strength  were  identified. 

The  manufacturing  and  analysis  work  done  to  date  has  demonstrated  the  feasi¬ 
bility  of  modeling  the  slip-cast  RBSN  process  with  the  Adaptive  Learning 
Network  methodology  and  is  viewed  as  the  first  iteration  in  the  optimization 
procedure  which  is  ultimately  aimed  at  finding  those  manufacturing  conditions 
which  will  produce  the  strongest,  most  consistent  material  strengths. 


2.  DATA  BASE 


Process  variables  and  resulting  part  strengths  were  recorded  for  35  different 
production  conditions.  Approximately  twenty  test  bars  were  manufactured  at 
each  condition,  and  the  bars  were  destuctively  stress-tested  to  determine 
their  strengths  in  terms  of  room  temperature  modulus-of-rupture  (MOR)  and 
strength  variance,  which  inversely  is  related  to  Weibul  Modulus.  For  the  work 
done  to  date,  the  strength  resulting  from  a  certain  production  condition  is 
taken  to  be  the  average  strength  of  the  twenty  test  bars.  A  listing  of  the 
data  base  is  presented  in  Appendix  1.  From  this  listing,  it  can  be  seen  which 
process  variables  were  measured  and/or  computed  as  well  as  the  numerical 
ranges  of  the  variables. 

The  raw  particle  size  distribution  (PSD)  data  and  the  sintering  and  nitriding 
temperature  histories  are  continuous  curves.  For  use  in  modeling,  discrete 
parameters  of  the  curves  must  be  computed.  It  is  as  crucial  to  compute 
"important"  parameters  from  continuous  data  as  it  is  to  instrument  any 
significant,  but  directly  measurable,  discrete  variable.  Discovering  useful 
waveform  parameters  is  often  complicated  due  to  the  very  large  number  of 
possibilities.  In  practice,  several  specific  parameters  from  each  of  several 
categories  are  computed  and  input  to  an  ALN  training  algorithm  which  selects 
the  most  useful  of  the  parameters  presented  to  it.  Based  on  which  variables 
are  selected,  knowledge  of  the  physical  process,  and  experience  from  other 
processes  which  may  have  some  similarities,  new  parameters  are  formulated  and 
tested  in  further  ALN  training.  Ultimately  it  is  desirable  to  reduce  any 
curve  to  three  or  four  key  descriptors. 


2.1  PARAMETERIZATION  OF  PARTICLE  SIZE  DISTRIBUTIONS 

To  begin  the  parameterization  of  the  PSD  curves,  the  cumulative  distribution 
curves  were  differentiated  and  scaled  to  provide  relative  particle  size 
density  functions. 

Twenty-six  parameters  (numbers  67  through  92  in  the  data  base)  were  computed 
from  the  PSD  curves.  Five  parameters  indicated  the  percentage,  by  weight,  of 
particles  greater  than  specified  sizes  (40,  20,  10,  5,  and  1  micrometer ( s )) . 

These  percentage  parameters  were  tantamount  to  evaluating  the  cumulative 
curves  at  specified  size  points.  Five  parameters  indicated  the  particle  size, 
in  logarithmic  form,  for  which  specified  percentage  levels  (20,  50,  80,  95, 
and  98)  were  achieved.  These  size  parameters  were  tantamount  to  evaluating 
the  cumulative  curves  at  specified  percentage  points.  The  next  parameter  was 
the  size  of  the  largest  particle  found  in  a  sample  of  the  powder.  Six  parame¬ 
ters  indicated  the  amount  of  weight  in  six  adjacent  size  bins  (.0  to  .3,  . 3  to 
1.0,  1.0  to  3.0,  3.0  to  10.0,  10.0  to  30.0,  and  above  30.0).  Five  parameters 
were  the  ratios  of  weights  in  the  various  bins.  The  final  four  parameters 
were  the  first  four  moment-generating  functions  indicating  average,  variance, 
skew,  and  kurtosis  of  the  density  function. 
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2. 1. 1  Modeling  Strength  as  a  Function  of  the  PSD  Parameters 


When  modeling  strength  from  the  above  original  PSD  parameters,  the  variables 
that  were  selected  by  the  netwrks  were  shape  parameters,  most  particularly 
the  second  and  fourth  moment  generating  functions  and  the  number  of  particles 
greater  than  40  pm.  The  second  moment  is  variance  which  indicates  the  spread 
or  width  of  the  distribution.  The  fourth  moment  is  kurtosis,  which  is  a  rough 
indicator  of  uni-  versus  bi-modality  of  a  curve. 

Investigation  into  the  models  shows  that  high  strengths  are  achieved  with  (a) 
a  broad  variance,  (b)  high  kurtosis,  i.e.,  uni-modality,  and  (c)  no  particles 
larger  than  40  ym.  These  analyses  are  borne  out  by  a  visual  inspection  of  the 
PSD  curves.  For  this  purpose,  the  PSD  curves  for  the  14  powders  are  shown  in 
Figure  2.  1.  The  curves  are  arranged  in  order  of  their  average  resulting 
strengths.  The  first  four  curves,  A^Bg— ^  (top  two),  A4Bg,  and  A2Bg  produced 
the  strongest  parts,  and  each  had  the  properties  selected  by  the  models.  The 
very  low  strength  curves  were  multimodal.  The  A2B7  curve  had  too  small  a 
spread. 


2.1.2  Hypothesis  of  an  Ideal  Particle  Size  Distribution 

From  the  above  results,  it  was  hypothesized  that  an  ideal  particle  size 
distribution  would  be  of  the  form: 

p(s)  =  As11  exp  ]  (1) 

where 

s  =  particle  size 

p(s)  =  particle  size  distribution 

A  =  normalization  scale  factor 

a  =  value  of  s  at  which  p(s)  peaks 

n  =  constant  controlling  the  rise  rate  of  p(s) 

m  =  constant  controlling  the  fall  rate  of  p(s) 

This  curve  is  unimodal,  which  accounts  for  the  kurtosis  constraint,  but  the 
coefficients  a,  n,  and  m  permit  a  range  of  shapes  to  fit  the  "good"  PSD 
curves. 


2.1.3  Fitting  the  Distributions 

The  first  step  in  testing  the  ideal  distribution  hypothesis  involved  fitting 
the  actual  particle-size  distributions  with  the  postulated  formula  and  finding 
the  coefficients  a,  n,  and  m  which  gave  the  least-squares  error  fit  to  each  of 
the  actual  distributions.  The  curve  fitting  approach  employed  a  gradient 
accelerated  search  to  find  the  optimum  values  of  the  coefficients  a,  n,  and  m. 


-I ^Letters  A,  B,  D,  F,  and  G  indicate  subprocesses:  A  -  starting  powder,  B  - 
powder  preparation,  D  -  slip  preparation,  F  -  sintering,  and  G  -  nitrid¬ 
ing.  The  subscripts  indicate  the  manufacturing  condition  number.  Thus 
A2Bg  indicates  the  second  starting  powder  processed  by  the  sixth  pro¬ 
cedure.  The  specific  condition  parameters  are  given  in  the  Appendix. 
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FIGURE  2.1:  ACTUAL  PARTICLE  SIZE  DISTRIBUTIONS  AND  CORRESPONDING 
VALUES  FOR  MEAN  STRENGTH,  SHOWN  IN  ORDER  OF  MEAN 
RESULTING  STRENGTH 
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A  gradient  search  was  selected  over  a  regression  fitting  approach  because  of 
the  difficulty  of  linearizing  the  theoretical  formula.  The  optimum 
coefficient  values  are  considered  to  be  those  which  produce  the  minimum  total 
k  deviation  between  the  postulated  and  actual  curves.  The  total  deviation,  or 

fitting  error,  is  defined  to  be  the  sum  of  the  squares  of  the  density 
differences  between  the  two  curves  at  each  point  along  the  size  axis. 

2.1.4  The  Gradient  Accelerated  Search  Algorithm 

The  starting  values  of  a,  n,  and  m  are  nominal  values  which  are  the  same  for 
each  actual  distribution.  This  starting  point,  and  the  resulting  error,  is 
initially  defined  to  be  the  best-to-date ,  and  the  error  gradient  is  computed 
at  that  point.  To  find  the  next  trial  point,  a  step  is  then  taken  from  the 
best-to-date  point  along  the  gradient  vector.  The  size  of  the  first  step  is 
user-specified  and  is  small.  The  resulting  theoretical  curve  is  compared  to 
the  actual  distribution,  and  the  fitting  error  is  computed.  If  the  resulting 
error  is  reduced,  the  trial  point  becomes  the  best-to-date. 

Rather  than  recomputing  the  gradient  and  moving  in  a  new  direction,  the  next 
step  is  taken  along  the  original  gradient.  The  step  is  taken  from  the  new 
best-to-date  point  but  the  step  size  is  doubled  (acceleration) .  Motion 
continues  along  the  same  gradient  line  with  the  step  size  continually  doubling 
until  the  performance  no  longer  improves.  The  search  then  returns  to  the  most 
recent  best-to-date  point,  computes  a  new  gradient  vector,  resets  the  step 
size  to  the  original  small  value,  and  proceeds  along  the  new  gradient  vector. 

The  search  is  stopped  when  the  first  (small)  step  along  a  newly  computed 
gradient  vector  does  not  produce  an  improvement.  The  search  is  then  within 
one  small  step  of  the  peak,  which  is  considered  to  be  sufficiently  close. 

The  acceleration  feature  of  the  search  allows  small  steps  to  be  taken  to 
pinpoint  the  peak,  but  avoids  the  excessive  computation  time  of  a  fixed-step- 
size  search  which  always  "creeps"  along. 

A  potential  problem  with  gradient  searches  is  that  they  can  get  "trapped"  on 
local  non-optimum  peaks.  It  has  been  established  that  the  search  space  for 
the  particle-size  distributions  is  unimodal  over  the  region  of  interest  so 
there  is  no  problem  of  identifying  the  global  peak. 

2.1.5  Fitting  Results 

i  The  theoretical  curves  fitted  to  each  of  the  fourteen  actual  particle-size 

distributions  are  shown  by  the  heavy  lines  in  Figure  2.2.  Table  1  gives  the 
values  of  the  coefficients  which  were  found,  the  residual  error  resulting  from 
the  optimum  fit,  and  the  strengths  of  the  test  bars  which  were  manufactured 
from  the  powder. 

An  analysis  of  these  data  shows  that  the  psd ' s  which  yield  high  strength 
typically  have  high  values  of  n  accompanied  by  low  values  of  m,  which 
indicates  that  excessive  amounts  of  powder  below  .3  yM  and  above  20.0  yM  are 
!  undesirable.  Also,  those  psd's  which  can  be  fitted  closely  by  the  theoretical 

curve  result  in  higher  ceramic  strengths  than  those  curves  which  cannot  be 
i  fitted  well.  This  tends  to  support  the  hypothesis  that  such  a  curve  is  ideal. 
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TABLE  2.1:  PARTICLE  SIZE  DISTRIBUTION 
FITTING  COEFFICIENTS 


Powder 

Number 

Condition 

Numbers 

Fitted 

a. 

Coefficient 

n ,  m. 

Values 

n/m 

RMS  Fitting 

Error 

Strengths  ( MOP. ) 
mean  min  max 

A1B1 

1-1 1 

3.63 

1. 09 

1.50 

0.73 

.  107 

26.5 

18.6 

33.  3 

A1B2 

12-15 

3.55 

0.90 

1.52 

0.59 

.103 

30.4 

26.8 

39.3 

A1B3 

16-17 

3.63 

1.  13 

1.47 

0.76 

.  124 

26.3 

24.6 

28.  0 

A1B4 

18-20 

3.89 

1.25 

1.52 

0.82 

.  113 

27.4 

25.9 

29.4 

A2B  1 

21 

5.88 

2.  36 

0.47 

5.02 

.  103 

33.  1 

33.  1 

33.  1 

A2B6 

22-25, 

30 

4.68 

3.49 

0.59 

5.91 

.089 

43.  1 

41.4 

45.0 

A2B17 

26 

6.61 

3.81 

3.  16 

1.21 

.091 

27.  7 

27.7 

27.7 

A3B 1 

27 

4.57 

3.05 

0.37 

0.68 

.068 

36. 1 

36. 1 

36.  1 

A3B5 

28-29 

8.  31 

0.  96 

3.  04 

0.31 

.  071 

27.6 

27.  1 

28.0 

a3B9 

31 

17.72 

.79 

1.80 

.44 

.015 

40.5 

26.0 

46.7 

A4B6 

32 

18.85 

.79 

1. 80 

.44 

•004 

46.  5 

40.5 

50.2 

A4B9 

33 

17.84 

3.91 

.75 

5.21 

.012 

36.3 

15.9 

56.0 

A4B6 

34 

17.  52 

5.  0 

.85 

5.88 

.012 

47.9 

39.0 

54.7 

A4B9 

35 

17.69 

5.0 

.67 

7.46 

.013 

44.0 

39.8 

48.5 

Fitted  Distribution:  p(s)  =  A  sn  exp [-n/m( s/a )m] 

s  =»  article  size 

p(s)  =  particle  size  distribution 

A  =  normalization  scale  factor 

a  =  value  of  s  at  which  p(s)  peaks 

n  =  constant  controlling  rise  rate  of  p(s) 

m  =  constant  controlling  fall  rate  of  p(s) 
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particle  size  process  condition 
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FIGURE  2.2:  ACTUAL  AND  FITTED  PARTICLE  SIZE  DISTRIBUTIONS  FOR  THE  14 
POWDERS  USED  IN  THE  35  EXPERIMENTS 


(continued) 
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PARTICLE  SIZE  PROCESS  CONDITION 
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FIGURE  2.2:  ACTUAL  AND  FITTED  PARTICLE  SIZE  DISTRIBUTIONS  FOR  THE  14 
POWDERS  USED  IN  THE  35  EXPERIMENTS 
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PARTICLE  SIZE  PROCESS  CONDITION 
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FIGURE  2.2:  ACTUAL  AND  FITTED  PARTICLE  SIZE  DISTRIBUTIONS  FOR  THE  14 
POWDERS  USED  IN  THE  35  EXPERIMENTS 


2. 2  PARAMETERIZATION  OF  THE  SINTERING  AND  NITRIDING  TEMPERATURE  PROFILES 


Seven  different  sintering  runs  and  three  different  nitriding  runs  were  used 
among  the  35  production  conditions.  The  temperature-versus-time  profiles  of 
the  sintering  and  nitriding  runs  are  shown  in  Figures  2.3  and  2.4  respec¬ 
tively,  with  the  resulting  average  strengths  are  shown  on  the  right  of  the 
Figures. 


2. 2. 1  Sintering  Profile  Parameters 

Fourteen  parameters  (numbers  113  through  136  in  the  data  base)  were  computed 
from  each  of  the  seven  sintering  temperature  profiles.  The  first  six  parame¬ 
ters  are  the  duration  times  that  the  sintering  temperatures  were  above  speci¬ 
fied  temperatures  (200°,  400°,  600°,  800°,  900°,  and  1000°C).  The  next  six 
parameters  are  the  degree  hours  above  the  (same)  specified  temperatures,  i.e., 
the  areas  under  the  curve  and  above  the  threshold  temperature.  The  final  two 
parameters  are  the  rise  and  fall  times  between  21°  and  900°  C. 


2.2.2  Nitriding  Profile  Parameters 

Fourteen  parameters  (numbers  141  through  154  in  the  data  base)  were  computed 
from  each  of  three  nitriding  temperature  profiles.  The  parameters  are  the 
same  as  the  sintering  variables  except  that  the  specified  temperature  levels 
were  400,  600,  800,  1000,  1200,  and  1300  °C.  No  rise  or  fall  times  wee 
computed  as  these  were  very  rapid  with  respect  to  the  overall  nitriding  times. 
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FIGURE  2.3:  TEMPERATURE  VS.  TIME  PROFILES  FOR  THE  SINTERING  PROCESS 
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FIGURE  2.3  (continued) 
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FIGURE  2.4:  TEMPERATURE  VS.  TIME 
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PROFILES  FOR  THE  NITRIDING  PROCESS 
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3.  ADAPTIVE  LEARNING  NETWORK  MODELING 


3. 1  ALN  MODELS 

Several  ALN  models  of  the  ceramics  manufacturing  process  were  generated  using 
the  Adaptronics,  Inc.  PNETTR(3)  model  synthesis  algorithm.  The  models 
discussed  here  were  trained  using  all  35  data  points  provided  by  AiResearch. 
Three  categories  of  models  were  developed:  (1)  strength  and  strength  variance 
were  modeled  as  a  function  of  the  independent  process  input  variables,  (2) 
strength  and  strength  variance  were  modeled  as  a  function  of  dependent 
intermediate  process  variables,  and  (3)  the  dependent  intermediate  process 
variables  were  modeled  as  a  function  of  the  independent  process  input 
variables. 

Figures  3.1a  through  3.12a  show  the  models  from  each  of  the  three  respective 
categories.  All  model  inouts  and  outputs  have  been  linearly  scaled  to  zero 
mean  and  unit  standard  deviation  (see  Appendix  1  for  scaling  factors)  to  allow 
an  evaluation  of  relative  variable  importance  by  comparison  of  coefficient 
magnitudes.  The  predominant  mathematical  terms  of  the  models,  in  unitized 
partial  derivative  form,  are  shown  along  with  the  network  block  diagrams. 
These  partial  derivatives  are  quantitative  estimates  of  the  relationships 
between  variables. 

Also  shown  in  each  ALN  figure  are  the  ranges,  R,  of  values  over  which  the 
models  were  trained,  the  standard  deviation,  S,  of  the  data,  and  the  RMS  error 
that  the  models  produced.  There  are  two  error  metrics.  The  first  is  the  RMS 
error,  e,  that  was  obtained  on  the  35  data  points  used  in  training,  and  the 
second  is  the  RMS  error,  E,  that  the  model  would  be  expected  to  make  on  new 
data  which  was  not  used  in  the  training  process,  E.  Model  usefulness  should 
be  judged  by  the  second  metric,  the  expected  error.  It  must  be  emphasized 
that  the  expected  error  is  a  valid  estimate  only  if  the  new  data  presented  to 
the  model  is  statistically  similar  to  the  original  training  data,  i.e.,  that 
the  values  of  the  input  and  output  variables  are  within  the  range  of  the 
training  data.  The  model  performance  measure,  P  =  1-E/S,  is  unity  minus  the 
ratio  the  expected  error  on  new  data  to  the  standard  deviation  of  the  original 
data.  If  this  number  is  equal  to  zero,  the  model  is  of  no  value?  if  it  is 
equal  to  unity,  the  model  is  a  perfect  predictor. 

To  provide  visual  insight  to  the  models,  contour  plots  of  the  ALN's  are 
presented  in  Figures  3.1b  through  3.12b.  These  diagrams  show  each  network's 
output  as  a  function  of  its  two  most  predominant  inputs.  If  a  model  has  more 
than  two  inputs,  the  values  of  those  inputs,  for  plotting  purposes,  are  held 
constant  at  their  mean  value.  The  curves  on  the  plots  show  contours  of  con¬ 
stant  model  output,  and  the  numbers  next  to  the  curves  indicate  the  model 
output  value. 

The  asterisks  on  the  contour  diagrams  show  the  locations  of  the  35  data  points 
used  in  the  model  synthesis.  Models  are  expected  to  be  most  accurate  in  the 
vicinity  of  the  data  and  less  accurate  further  away  from  the  data.  There  are 
not  always  35  asterisks  on  each  plot.  In  many  cases  several  observations  had 
identical  values  for  the  two  parameters  being  plotted,  so  an  asterisk  may 
represent  several  points. 
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3.2  MODEL  INTERPRETATIONS 


Investigation  of  the  model  structures  leads  to  the  following  hypotheses  about 
the  slip-cast,  reaction-bonded  silicon-nitride  manufacturing  process.  Because 
of  the  small  amount  of  data  used  in  the  model  synthesis,  these  interpretations 
should  be  viewed  as  no  more  than  hypotheses  about  the  process.  It  is  intended 
that  the  interpretations  be  used  only  as  guides  for  further  experimentation 
and  not  as  definitive  statements  about  the  chemical  process. 

3. 2.  1  Mean  Strength  Modeled  as  a  Function  of  the  Independent  Variables; 
(Figures  3.1  and  3.2) 

The  presence  of  oxygen  in  the  starting  powder  has  a  detrimental  effect  on 
strength;  therefore  a  low  amount  (less  than  0.5  percent)  of  oxygen  in  the 
starting  powder  is  desirable. 

The  use  of  larger  amounts  of  media  quantity  (above  10  Kg  Al2C>3)  in  the  powder 
preparation  has  a  positive  effect  on  the  ultimate  ceramic  strength. 

Increasing  the  number  of  particles  greater  than  40  m  generally  reduces 
strength.  There  are  two  notable  exceptions  to  this  trend  as  shown  by  the  A2B1 
and  A3B1  powders  in  Figure  3.2b.  Both  of  these  had  approximately  one  percent 
of  the  particles  greater  than  40  m  yet  still  achieved  actual  strengths  of  33 
and  36  ksi  respectively,  and  these  two  data  points  account  for  the  cruadratic 
term  in  model.  But  the  majority  of  the  data  lies  to  the  left  on  the  plot,  with 
less  than  0.33%  of  particles  greater  than  40  m,  and  in  this  region  it  is 
apparent  that  smaller  percentages  of  large  particles  contribute  to  higher 
strength. 

Increasing  the  standard  deviation,  or  spread,  of  the  particle  size  distri¬ 
bution  has  a  minor  positive  effect  on  strength.  Very  narrow  distributions 
should  be  avoided. 


3. 2. 2  Strength  Variance  Modeled  as  a  Function  of  the  Independent  Variables; 
(Figure  3.3) 

Decreasing  the  rise  coefficient  (n)  of  the  fitted  particle  size  distribution 
tends  to  decrease  strength  variance.  Lowering  n  corresponds  to  increasing  the 
proportion  of  smaller  particles  in  the  overall  size  distribution  and  is  in 
keeping  with  the  requirement  for  a  broad  particle  size  distribution. 

Increasing  the  coefficient  of  skewness  of  the  particle  size  distribution  tends 
to  reduce  the  strength  variance.  Thus,  though  the  total  distribution  should 
be  relatively  broad  and  its  mean  shifted  twoard  the  smaller  sizes,  its  shape 
should  be  skewed  toward  the  larger  sizes. 

Shorter  sintering  times  (less  than  10  hours)  at  temperatures  greater  than 
900°C  appear  to  decrease  the  strength  variance. 
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(R)  Range  of  the  Data  :  29.1  (ksi)  Min  =  18.69, 

Max  =  47.80 

(S)  Standard  Deviation  :  7.82  (ksi) 

(e)  RMS  Error  on  Training  Data  Base  s  3.88  (ksi) 

(E)  Expected  RMS  Error  on  New  Data  :  4.6  (ksi) 

(P)  Model  Performance  Measure  ( 1-E/S)  :  0.41 


Partial  Derivatives  of  Mean  MOP  with  Respect  to  the  Model  Input  Variables: 

Oxygen  in  Starting  Powder  (%)  :-7.71  (ksi/%) 

Media  Quantity  (Al^O^/kg)  :  1.03  (ksi/kg) 


FIGURE  3.1a:  ALN  MODEL  PREDICTING  STRENGTH  AS  A  FUNCTION  OF  INDEPENDENT 

VARIABLES 
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(P) 
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( 1 -E/S ) 

:  .34 

Min 
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18.69 

47.80 


Partial  Derivative  of  MOR  Mean  with  Respect  to  the  Model  Input  Variables: 

Percent  Greater  than  40  ym  :-37.90  (ksi/%) 

Standard  Deviation  of  log  of  PSD  :  3.15  (ksi/  M) 


FIGURE  3.2a:  ALN  MODEL  PREDICTING  STRENGTH  AS  A  FUNCTION  OF  PARTICLE  SIZE 
DISTRIBUTION  PARAMETERS 
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(R)  Range  of  the  Data 

(S)  Standard  Deviation 

(e)  RMS  Error  on  Training  Data  Base 
(E)  Expected  RMS  Error  on  New  Data 
(p)  Model  Performance  Measure  (1-E/S) 


10.82  (ksi);  Min  =  2.84 

Max  =  13. 66 

2.19  (ksi) 

1.22  (ksi) 

1.70  (ksi) 

0.  22 


Partial  Derivatives  of  MOR  Standard  Deviation  with  Respect  to  the  Model  Input 
Variables: 

Fitted  PSD  Rise  Coefficient  ( B)  :  .81  (ksi) 

Coefficient  of  Skewness  :-.73  (ksi) 

Sintering  Time  Greater  Than  900«C  Hrs  :  .68  (ksi/°C  Hrs) 


ALN  MODEL  PREDICTING  STRENGTH  VARIANCE  AS  A  FUNCTION  OF  INDE 
PENDENT  VARIABLES 


FIGURE  3.3a: 


Mean  Strength  Modeled  as  a  Function  of  the  Intermediate  Process 
Variables:  (Figure  3.4) 


High  strength  results  from  high  nitrided  density  and  high  weight  gain  during 
the  nitriding  process.  As  would  be  expected  from  a  chemical  analysis,  the 
weight  gain  should  approach  the  theoretical  maximum  of  approximately  62%.  The 
nitrided  density  would  ideally  be  above  2.75  gm/cm^. 


Low  percentages  of  Alpha-Silicon-Nitride  (less  than  75  percent)  in  the  final 
analysis  yield  high  strength. 


3.2.4  Strength  Variance  Modeled  as  a  Function  of  the  Intermediate  Process 
Variables  (Figure  3.5) 

Decreasing  the  ratio  of  Silicon  Oxy-Nitride  to  Alpha-Silicon-Nitride  tends  to 
decrease  strength  variance. 

Decreasing  the  ratio  of  Beta-  to  Alpha-Silicon  Nitride  appears  to  decrease 
strength  variance  if  the  ratio  is  .25  or  higher  to  start  with.  Decreasing 
this  ratio  implies  increasing  Alpha  which  (from  Section  3.2.3)  would  reduce 
strength,  so  it  appears  that  there  is  a  small  tradeoff  between  strength  and 
variance.  Since  low  Alpha  has  a  greater  positive  effect  on  strength  than  it 
has  an  adverse  effect  on  variance,  it  is  recommended  that  Alpha  be  minimized. 

Decreasing  the  1/30  viscosity  appears  to  reduce  strength  variance  to  a  small 
degree. 


3.2.5  Intermediate  Process  Variables  as  a  Function  of  the  Independent 
Variables: 

Slip  pH:  (Figure  3.6) 

Slip  pH  appears  to  be  highly  nonlinearly  dependent  upon  percent  of  solids 
in  the  slip  and  the  percent  of  manganese  in  the  starting  powder.  The 
model  is  quite  biased,  however,  by  observation  number  26,  which  came  from 
the  A2B.17  powder,  had  an  extremely  low  solids  content  of  52%,  and  resul¬ 
ted  in  a  moderately  low  mean  strength  of  27.7  ksi. 

Green  Density:  (Figure  3.7) 

Increasing  the  percentage  of  slip  additive  NH^OH  results  in  increased 
Green  Density. 

For  values  between  0.00  and  0.04  percent,  the  amount  of  deflocculant  does 
not  significantly  affect  green  density,  but  above  about  0.04  percent, 
additional  defloculant  appears  to  increase  Green  Density. 

Increasing  the  quantity  of  milling  medium  ALjOj  and  increasing  the  slip 
aging  time  both  cause  minor  increases  in  the  Green  Density  of  the  parts. 
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(R)  Range  of  the  Data 

(S)  Standard  Deviation 

(e)  RMS  Error  on  Training  Data  Base 
(E)  Expected  RMS  Error  on  New  Data 
(P)  Model  Performance  Measure  (1-E/S) 


29.1  (ksi)  Min  =  18.69, 
Max  =  47.80 

7.82  (ksi) 

3.59  (ksi) 

4.72  (ksi) 

0.  40 


Partial  Derivatives  of  Mean  Strength 

Nitrided  Density  (gm/cm3) 

Weight  Gain  (%) 

Alpha  (Rel.  %) 


with  Respect  to  Model  Input  Variables: 

:  24.4  (ksi/ (gm/cm3) 

:  1.63  (ksi/%) 

0.80  (ksi/%) 


FIGURE  3.4a:  ALN  MODEL  PREDICTING  MEAN  STRENGTH  AS  A  FUNCTION  OF  INTER¬ 
MEDIATE  PROCESS  VARIABLES 
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Viscos  1/30 
Beta/Alpha 


Si 2ON2/Alpha 


(R)  Range  of  the  Data 

( S )  Standard  Deviation 

(e)  RMS  Error  on  Training  Data  Base 
(F)  Eexpected  RMS  Error  on  New  Data 
(P)  Model  Performance  Measure  (1-E/S) 


10.82  (ksi) 

2.19  (ksi) 
1.23  (ksi) 
1.63  (ksi) 
0.26 


Min  =  2.84 

Max  =  13.66 


Partial  Derivatives  of  MOR  Standard  Deviation  with  Respect  to  the  Model  Input 
Variables: 


Ratio  Si2ON2/Alpha  (%/%) 
Ratio  Beta/Alpha  (%/%) 
Viscosity  1/30  (CPS) 


33.95  (ksi) 

8.76  (ksi) 

0.02  (ksi/ CPS ) 


FIGURE  3.5a:  ALN  MODEL  PREDICTING  STRENGTH  VARIANCE  AS  A  FUNCTION  OF  INTER¬ 
MEDIATE  PROCESS  VARIABLES 
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x5- 

[  2. 17  *  2. 58x^  “  1 • 52xgg  ^  5*22x^Xgg 
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-  I.I8X5  -  1.31Xgg 

X1 05 
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(R) 

Range  of  the  Data 

2.2 

(-) 

(S) 

Standard  Deviation 

0.545 

(-) 

(e) 

RMS  Error  on  Training 

Data 

0.  19 

(-) 

(E) 

Expected  Error  on  New 

Data 

0.34 

(-) 

(P) 

Model  Performance  Measure  (1-E/S) 

.37 

=  4. 9,  Max  ■  5.9 


Partial  Derivatives  of  Ph  with  Respect  to  the  Model  Input  Variables: 


%  Manganese  in  Starting  Powder  (%) 
Solids  Content  in  Slip  (%) 


:-69.66  (-/%) 
.21  (-/%) 


FIGURE  3.6a:  ALN  MODEL  PREDICTING  SLIP  pH  AS  A  FUNCTION  OF  INDEPENDENT 

VARIABLES 
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Max  =  1 . 
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:  0. 40 

Partial  Derivatives  of 

Green 

Density 

with 

Respect 

to  the 

Variables: 


Additive  NH4OH  (%  Wght) 
Deflocculant  (%  Wght) 
Media  Quantity  (kg) 

Slip  Aging  Time  (Days) 


2.50  ( gm/cm3 /%) 

2.25  (gm/cm3/%) 

0.01  (gm/cm3/Kg) 
0.003  (gm/cm3/day) 


FIGURE  3.7a:  ALN  MODEL  PREDICTING  GREEN  DENSITY  AS  A  FUNCTION 
VARIABLES 
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Weight  Gain:  (Figure  3.8) 


Weight  gain  increases  proportionally  to  the  particle  size  distribution 
ratio  [(Weight  Bin3  plus  Weight  Bin5)  divided  by  Weight  Bin4].  This 
ratio  reflects  the  width  of  the  PSD,  indicating  that  a  broad  PSD  yields 
high  weight  gain. 

Increasina  the  percentage  of  iron  in  the  starting  powder  tends  to 
decrease  the  weight  gain. 

Decreasing  the  slip  aging  time  appears  to  cause  a  minor  increase  in 
weight  gain. 

Nitrided  Density;  (Figure  3.9) 

From  the  quadratic  nature  of  the  model,  it  appears  that  weight  gain  is 
relatively  insensitive  to  the  amount  of  deflocculant  if  the  deflocculant 
is  less  than  about  0.04%  of  the  slip,  but  above  0.04%,  additional 
deflocculant  increases  weight  gain.  But  because  of  the  very  limited  data 
above  0.03%,  this  conclusion  must  be  considered  to  be  very  weak. 

Decreasing  the  skewness,  i.e.,  increasing  the  amount  of  small  particle 
sizes,  of  the  particle  size  distribution  tends  to  increase  the  nitrided 
density. 

Increasing  the  solids  content  of  the  slip  appears  to  cause  a  minor 
increase  in  the  nitrided  density. 

Alpha :  (Figure  3.10) 

Decreasing  the  amount  of  FejO-j  additive  in  the  slip  decreases  the  propor¬ 
tion  of  Alpha-Silicon-Nitride  in  the  final  analysis. 

Decreasing  the  amount  of  oxygen  in  the  starting  powder  decreases  the 
proportion  of  Alpha-Silicon  Nitride. 

Decreasing  the  slip  temperature  causes  a  minor  decrease  in  the  percentage 
of  Alpha-Silicon-Nitride  in  the  final  analysis. 

Beta:  (Figure  3.11) 

A  lower  percentage  of  Manganese  in  the  starting  powder  results  in  a 
higher  proportion  of  Beta-Silicon-Nitride  in  the  final  analysis. 

Lower  percentages  of  deflocculant  in  the  slip  preparation  result  in 
higher  proportions  of  Beta-Silicon-Nitride. 

Decreasing  the  20th  percentile  size  of  the  PSD  (i.e.,  increasing  the 
number  of  small  particles)  tends  to  increase  slightly  the  proportion  of 
Beta-Silicon-Nitride  in  the  final  product. 

Decreasing  the  slip  temperature  has  a  minor  tendency  to  increase  the 
final  percentage  of  Beta-Silicon-Nitride. 


X155  Gain 
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Range  of  the  Data 
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7.7  (%)  Min  =  53.5,  Max  =  61.2 

1.78  (%) 

0.89  (%) 

1. 19  (%) 

0.33 


Partial  Derivatives  of  %  Weight  Gain  with  Respect  to  the  Model  Input  Variables 


Ratio  ( Bin3  +  Bin5)/Bin4 
Fe  in  Starting  Powder  ( % ) 
Slip  Aging  Time  (days) 


4.01  (%/%) 
-1.71  (%/%) 
-0.14  (%/day) 


FIGURE  3.8a:  ALN  MODEL  PREDICTING  WEIGHT  GAIN  AS  A  FUNCTION  OF  INDEPENDENT 
VARIABLES 
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Partial  Derivatives  of  Nitrided  Density  with  Respect  to  the  Model  Input 
Variables 


Deflocculant  (%  of  Wght) 
Coefficient  of  Skewness 
Solids  Content  (%  of  Weight) 


2.64  ( gm/ cmJ/% ) 
-0.09  (gm/cm3/-) 
0.05  (gm/cm3/%) 


FIGURE  3.9a:  ALN  MODEL  PREDICTING  NITRIDED  DENSITY  AS  A  FUNCTION  OF  INDE¬ 
PENDENT  VARIABLES 
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FIGURE  3.9b:  CONTOURS  OF  NITRIDED  DENSITY  PLOTTED  AS  A  FUNCTION 
OF  INDEPENDENT  VARIABLES 
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Partial  Derivatives  of  Alpha  with  Respect  to  the  Model  Input  Variables 

Slip  Additive  Fe203  (%)  :  3.72.0  {%/%) 

Oxygen  Content  in  Starting  Pwdr.(%)  :  3.78  (%/%) 

Slip  Temperature  (°F)  :  0.84  (%/%) 


FIGURE  3.10a:  ALN  MODEL  PREDICTING  ALPHA  AS  A  FUNCTION  OF  INDEPENDENT 
VARIABLES 
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PERCENT  OXYGEN  IN  STARTING  POWDER 


(R)  Range  of  the  Data 

(S)  Standard  Deviation 

(e)  RMS  Error  on  Training  Data  Base 

(E)  Expected  Error  on  New  Data 

(P)  Model  Performance  Measure  (1-E/S) 


Partial  Derivatives  of  Beta  with  Respect 

%  Manganese  in  Starting  Powder 
Deflocculant  (%  of  Weight) 

20  Percentile  Size  (Log) 

Slip  Temperature  (°F) 


21.4  (Rel.  %)  Min  =  7. 8, Max  =  29.2 

6.1  (Rel.  %) 

2.78  (Rel.  %) 

3.91  (Rel.  %) 

0.36 


to  the  Model  Input  Variables: 

:  -347.7  (%/%) 

:  -262.  3  (%/%) 

:  -  1.39  (%/-) 

:  -  0.98  (%/°F) 


FIGURE  3.11a:  ALN  MODEL  PREDICTING  BETA  AS  A  FUNCTION  OF  INDEPENDENT 
VARIABLES 
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Silicon-Oxy-Nitride  (Si2ON2):  (Figure  3.12) 


Titanium  in  the  starting  powder  increases  the  percentage  of  silicon-oxy¬ 
nitride  in  final  material. 

Adding  Fe202,  however,  appears  to  reduce  slightly  the  amount  of  silicon- 
oxy-nitride. 

3.3  SUMMARY  OF  KEY  EFFECTS 

The  key  hypotheses  obtained  from  the  networks  are  diagrammed  in  Figures  3. 13a 
and  b.  High  average  strength  results  primarily  from  high  weight  gain  and  high 
nitrided  density,  which  in  turn  are  achieved  with  a  broad,  unimodal  particle 
size  distribution  with  with  a  small  mean  size  and  with  low  iron  in  the  start¬ 
ing  powder.  low  alpha-silicon-nitride  also  improves  strength  and  is  achieved 
by  minimizing  oxygen  in  the  starting  powder  and  by  using  smaller  amounts  of 
ferris  oxide  additive  in  the  slip  preparation. 

Low  strength  variance,  i.e.,  high  Weibull  modulus,  is  achieved  by  minimizing 
silicon-oxy-nitride  and  to  a  lesser  extent  by  using  smaller  particle  sizes. 
The  silicon-oxy-nitride  appears  to  be  increased  most  by  titanium  in  the 
starting  powder. 

Production  conditions  yielding  high  average  strengths  also  yield  fairly  con¬ 
sistent  strengths,  i.e.,  high  Weibul  modulus,  while  low  average  strengths  are 
generally  accompanied  by  large  variation  in  part  strengths. 

Within  the  range  of  the  35  data  observations  used  in  this  analysis,  it  does 
not  appear  that  slip  pH,  green  density,  or  sintering  parameters  have  a 
significant  impact  on  strength  or  strength  variance.  Though  nitriding  param¬ 
eters,  such  as  temperatures,  times  or  nitrogen  pressures,  were  not  selected  by 
the  networks,  it  must  be  pointed  out  that  there  was  not  significant  variation 
of  those  parameters  in  the  data  base.  Further  data  should  be  collected  to 
determine  the  affects  of  nitriding  on  strength  and  strength  variance. 
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%Ti  Xg 

Fe2°3  x101 
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(R)  Rancje  of  the  Data 

(S)  Standard  Deviation 

(e)  RMS  Error  on  Training  Data  Base 
(E)  Expected  RMS  Error  on  New  Data 
(P)  Model  Performance  Measure  Cl -E/S) 


5.5  (%)  Min  =  0.0,  Max  =  5.5 

1.29  (%) 

1.0  (%) 

1.06  (%) 

0.  18 


Partial  Derivatives  OF  %  Si2ON2  with  Respect  to  the  Model  Input  Variables 

%  Titanium  in  Starting  Powder  :  56.12  (%/%) 

%  Slip  Additive  (Fe203)  :  -41.28  (%/%) 


FIGURE  3.12a:  ALN  MODEL  PREDICTING  %  Si2ON2  AS  A  FUNCTION  OF  INDEPENDENT 
VARIABLES 


I 


I 


3-28 


Major  Positive  Effect 
Minor  Positive  Effect 
Minor  Adverse  Effect 

FIGURE  3.13a:  FLOW  CHART  OF  PREDOMINANT  SLIP  CAST  RSSN  PROCESS  EFFECTS 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  Adaptive  Learning  Network  methodology  has  been  demonstrated  to  be  a 
powerful  tool  for  the  analysis  of  a  slip-cast,  reaction-bonded,  silicon- 
nitride  manufacturing  process.  Though  the  data  is  limited  (only  35  manu¬ 
facturing  variations  were  available)  ,  many  trends  have  been  identified  which 
will  be  useful  in  guiding  a  continued  search  for  the  optimum  manufacturing 
conditions. 

The  models  developed  to  date  do  not  exhibit  clear  cut  peaks  showing  optimum 
values  parameter  settings;  rather,  they  show  trends  which  suggest  that  further 
variations  of  the  parameters  will  yield  improved  material  properties.  It 
appears  that  strengths  well  above  48  ksi  (the  strongest  achieved  to  date)  are 
possible. 

There  was  little  variation  of  the  nitriding  parameters  in  the  given  data  base, 
so  the  impacts  of  nitriding  on  material  strength  could  not  be  estimated.  In 
future  work  it  is  recommended  that  data  be  collected  for  a  wide  range  of 
nitriding  conditions. 

It  is  also  recommended  that  future  work  address  high  temperature  strengths  as 
well  as  room  temperature  strengths,  so  that  strengths  may  be  optimized  for  the 
operational  environment  of  the  RBSN  materials 
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APPENDIX  1 


SLIP-CAST  REACTION-BONDED  SILICON-NITRIDE  DATA  BASE 

Notes : 

(1)  The  data  base  is  in  two  parts.  Section  1  presents  the  following 

statistics  for  each  variable  in  the  data  base:  mean,  standard 

deviation,  range,  minimum  value,  and  maximum  value.  Section  2 
presents  the  35  individual  values  of  each  variable. 

(2)  In  Section  2,  an  asterisk  next  to  a  zero  entry  indicates  that  no 
measurement  was  made  for  that  entry. 

(3)  Particle  size  data  is  given  in  the  following  log  form: 

s  =  size  in  micrometers 
x  =  log  size 
x  =  10  log  (10s) 
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